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ABSTRACT 



The multiple protostellar system L1551 IRS5 exhibits a large-scale bipolar 
molecular outflow that spans ~1.5 pc on both the NE (redshifted) and SW 
(blueshifted) sides of the system. We have studied this outflow within ~4000 AU 
of its driving source(s) with the SubMillimeter Array. Our C0(2-l) image at 
~4" (~560 AU) resolution reveals three distinct components: 1) an X-shaped 
structure spanning ~20" from center with a similar symmetry axis and velocity 
pattern as the large-scale outflow; 2) an S-shaped structure spanning ~10" from 
center also with a similar symmetry axis but opposite velocity pattern to the 
large-scale outflow; and 3) a compact central component spanning ~1.4" from 
center again with a similar symmetry axis and velocity pattern as the large-scale 
outflow. The X-shaped component likely comprises the limb-brightened walls 
of a cone-shaped cavity excavated by the outflows from the two main proto- 
stellar components. The compact central component likely comprises material 
within this cavity newly entrained by one or both outflows from the two main 
protostellar components. The S-shaped component mostly likely comprises a 
precessing outflow with its symmetry axis inclined in the opposite sense to the 
plane of the sky than the other two components, taking the S-shaped compo- 
nent out of the cone-shaped cavity along most if not all of its entire length. 
This outflow may be driven by a recently reported candidate third protostel- 
lar component in L1551 IRS5, whose circumstellar disk is misaligned relative to 
the two main protostellar components. Gravitational interactions between this 
protostellar component and its likely more massive northern (and perhaps also 
southern) neighbor (s) may be causing the circumstellar disk and hence outflow 
of this component to precess. 
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Introduction 



Most near-solar-mass stars are 



( iDuquennoy fc Mayor 



1991 



3orn a s members of binary or multiple systems 



Mathieu 



19941 ). The formation of such systems, although the 
preferred mode of star formation at least for near-solar-mass stars, is poorly understood. 
The main competing theoretical models are fragmentation of dense molecular condensations 
to form systems comprising multiple protostars, or capture of originally single protostars 
still embedde d in their individual condensations to form multiple systems (e.g., review by 



Tohline 



20021 ) 



The most direct way of distinguishing between these two possibilities is to compare 
for a given system the properties of its individual protostars (i.e., orientation of their 
circumstellar disks and their orbital motion) with the properties of its surrounding 
molecular condensation (i.e., orientation and rotation). Such a study of the low-mass 
protostellar system L1551 IRS5 provides the most direct evidence yet that its two main 
protostellar co mponents formed via th e fragmentation of their surrounding parent molecular 



condensation (ILim fc Takakuwa 



20061 ). For many systems (especially those with closely 



separated components), however, the properties of their individual protostellar components 
are difficult to measure until the Extended Very Large Array (EVLA) or Atacama Large 
Millimeter and subMillimeter Array (ALMA) is completed. Indeed, at the present time, 
even the number of protostellar components present in a given system is not always clear. 

Bipolar molecular outflows, found ubiquitously around protostars, are more easily 
accessible, and can be used to indirectly probe the properties of multiple protostellar 
systems. For example, such outflows can reveal the presence of currently undetectable or 
spatially unresolved protostellar components through their individual outflows. Outflows 
that wind in space with a helical pattern may have driving sources whose circumstellar disks 
and therefore jets are precessing due to gravitational interactions with closely-separated 
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protostellar companions. The degree of collimation and/or dynamical age of the outflow, 
as well as the opening angle of the outflow cavity at its base, can provide a guide to the 
evolutionary stage of the driving source. 

Here, we study the bipolar molecular outflow associated with L1551 IRS5 at a high 
angular re solution with an in terferometric array. This system was first found as an infrared 



source by 



Strom et al. 



(Il976l ) in the nearby (distance ~140 pc) L1551 molecular cloud 
located in the constellation Taurus. L1551 IRS5 was the first protostellar source from 
which a bipolar molecular outfl ow, detected in si ngle-dish observations of the CO (1-0) and 



CO (2-1) lines, was recognized (ISnell et al. 



19801 ). Subsequent single-dish observations in 



CO (1-0), CO (2-1), and CO (3-2) at higher angular resolutions have revealed the overall 
morphology and kinema tics of this and a number of other outflow s in the L1551 star-formin g 



region on parsec scales fjUchida et al. 



987 



Bachiller et al 



19941 



Moriartv-Schieven et al. 



1987 



Moriartv-Schieven fc Snell 



1988 



Moriartv-Schieven et al 



2006 



Stojimirovic et al. 



20061 ). The most prominent of these outflows is centered on L1551 IRS5, and extends ~1.5 
pc towards the NE (redshifted lobe) and SW (blueshifted lobe) directions at a position angle 
of ~50°. Each outflow lobe has a U-shaped structure, together with knotty components 
having distinct velocities. The other outflows detected in the L1551 cloud have been 
associated with other driving sources except for a collimated redshifted outflow with one 
end near L1551 IRS5 and extending ~2 pc to the west, the origin of which is not known. 

Apart from the large-scale bipolar molecular outflow, L1551 IRS5 is surrounded by a 
relatively c ompact molecu l ar con densation (radius of ~1200 AU in C^^O) that has been best 



studied by 



Momose et al. 



(119981 ). This condensation is elongated in the NW-SE direction, 
perpendicular to the major axis of the large-scale bipolar molecular outflow. The measured 
kinematics indicate a flat tened envelope that is rot ating and contract ing, as also suggested 



in earlier observations by 



Ohashi et al. 



fll996h and 



Saito et al. 



( 119961 ). This condensation 
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has been observed in C S (7-6), which t r aces •w armer and denser gas closer to the central 



protostellar system, by 



Takakuwa et al. 



(|2004J ). In the inner region of the envelope, the 



kinematics is dominated by rotation. 

The first ind i cation that L1551 IRS5 may be a binary system was reported by 



Sieging fc Cohen 



(119851 ). who showed that this system comprised two sources separated by 



■^0.3" (~42 AU) in the north-south direction at a wavelength of 2 cm. Subsequent studies 



resolved these two sources into twiii ionize d jets closely aligned with the larg e-scale bipolar 

s3 



molecular outfiow fiRodrfguez et al. 



2003bD. 



Looney. Mundy fc WelchI (119971 ) showed from 



high angular-resolution (~0'.'5) observations in the 2.7-mm continuum that L1551 IRS5 
contains two circur nstellar dust disk s surro unded by a circumbinary dust disk and/or 



extended envelope. 



Rodriguez et al. 



( 1l998l ) spatially resolved the two circumstellar dust 



disks (along their major axes) for the first t ime at 7 mm, and showed that they are centered 



on the above-mentioned twin ionized jets. 



Rodriguez et al. 



(l2003al ) showed that these 



two protostellar compon ents were (likely) orbiting each other in a clockwise direction. 
Lim fc Takakuwa (i2006i) spatially resolved the two circumstellar disks along both their 
major and minor axes, and showed that these disks are aligned with each other as well 
their surrounding rotating and fiattened molecular condensation (pseudodisk). They also 
showed that the two components were orbiting each other in the same direction as the 
rotation of the pseudodisk. All these properties are consistent with the notion that the 
two protostellar components formed as a result of fragna e ntatio n in the central region of 



their surrounding pseudodisk. Finally, 



Lim fc Takakuwal ( 120061 ) found a candidate third 



protostellar component located (in projection) near the northern protostellar component. 
The circumstellar disk of this third component is misaligned relative to the circumstellar 
disks of the two main components as well as their surrounding pseudodisk. 



We report in this paper high angular resolution (~4") observations in C0(2-l) 
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centered on the base of the large-scale bipolar molecular outflow from L1551 IRS5. There 



is, to the best of our knowledge, only one other published inter 



molecular outflow in CO, that presented by 



Barsony et al. 



'erometric image of this 



fll993h in CO (1-0) who studied 



relatively extended (> 1') features in the blue outflow lobe. In our study, we investigate the 
structure and kinematics of the outflow components detected within ~30" of L1551 IRS5, 
their connection with the large-scale bipolar molecular outflow imaged with single-dish 
telescopes, and the likely nature of their driving sources. 

Readers interested in how the observations were conducted and data reduced should 
now proceed to §2. Those interested only in the results of the C0(2-l) and simultaneous 
1.3 mm continuum measurements can skip ahead to §3. In §4, we show that our C0(2-l) 
maps likely traces three separate outflow components: the limb-brightened walls of an 
outflow cavity, material in the cavity entrained by one or both jets from the two main 
protostellar components, and a candidate third precessing bipolar outflow. Readers 
interested only in a concise summary of our results and interpretation can proceed directly 
to 55. 



2. Observations and Data Reduction 



Observations of L1551 IRS5 in C0(2-l) and 1.3 mm continuum were carried out usine 

n 

the SMA|^ in its compact configuration on 2003 December 7. A description of the SMA 



can be found in 



Ho et al. 



(120041 ). The parameters of the observations are summarized in 



Table 1. Seven of the eight antennas of the SMA were available for the observations. We 



^The Submillimeter Array (SMA) is a joint project between the Smithsonian Astrophys- 
ical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics and is 
funded by the Smithsonian Institution and the Academia Sinica. 
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used the stronger quasar 0423-013, lying 19.6° from L1551 IRS5, for amplitude calibration, 
and the weaker quasar 3C 120, lying 12.8° from L1551 IRS5, for phase calibration. The flux 
density of 0423-013 was 9.4 Jy and that of 3C 120 was 3.3 Jy as measured with respect to 
Uranus, which served as the absolute flux calibrator. 

The SMA has a double-sideband receiver in the 230 GHz (1.3 mm) band with a 
bandwidth of 2 GHz in each sideband. The lower sideband was tuned to a central frequency 
of 231.3 GHz, and the upper sideband 241.3 GHz. The signal in each sideband is fed 
into a correlator that distributes this signal to twenty- four spectral windows ("chunks" in 
the SMA nomenclature). Six of the chunks were divided into 512 channels, which in the 
CO (2-1) line results in a velocity resolution of 0.265 km s~^ over a velocity range of 135 
km s~^. The remaining chunks were divided into 128 channels, and following calibration 
vector-averaged to make a single continuum channel at 1.3 mm with a total bandwidth of 
4.0 GHz. The minimum projected baseline in our observation was ~8 kA, allowing us to 



[e.g., see 



^30% of the flux 


Tom 


Wilner & Welch 


1994) 



We calibrated the raw visibility data using MIR, which is an IDL-based data reduction 
package adopted for the S MA from the MMA software package originally developed for 



OVRO f Scoville et al 



19931 ). We adopted antenna-based calibration, which provides the best 
sensitivity, as there were no significant baseline-based errors in our data. The calibrated 
visibility data were Fourier-transformed to produce DIRTY images, and the point spread 
function of the te lescope deconvolved (CLEANed) from these images using MIRIAD 



( ISault et al. 



19951 ) to produce the final maps. A ROBUST parameter of 0.5 was adopted. 



which provided the best compromise between sensitivity and angular resolution. 

To make the CO (2-1) channel maps, we first subtracted the continuum emission from 
the visibility data as derived from a fit to the line-free channels. Unlike in the continuum 
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where the emission is centrally concentrated, in CO (2-1) the emission in many channels is 
distributed in a complex manner over an extended region resulting in complicated sidelobe 
patterns. To CLEAN the channel maps, we first selected (i.e., place boxes enclosing) 
features that we believed to be real. Through trial and error, we converged on features that 
changed smoothly in structure between neighboring channels, as would be expected for real 
features. An examination of the residual maps (i.e., after the candidate real features and 
their sidelobes had been removed) revealed no systematic pattern (i.e., only random noise), 
indicating that all the real feature had indeed been correctly recovered. Thus, we believe 
that any artificial components incorrectly selected for CLEANing are much weaker than the 
real components in our maps. 



Results 



3.1. 1.3 mm Continuum 



Figure 1 shows our 1.3-mm continuum map of L1551 IRS5 at an angular resolution 
of 4'.'66 X 2'.'33. The continuum emission peaks at the location of the protostellar systera , 



as se en in images a t comparable angular resolution at 3 mm (ILooney. Mundy fc Welch 



19971 ) and 0.8 mm (iTakakuwa et al. 



20041 ). We fitted a two-dimensional Gaussian to the 
continuum source to derive a total fiux density of 1.20±0.01 Jy and a deconvolved size at 
full-width half-maximum (FWHM) of (I'.'S ± O'.'l) x (O'.'S ± 0'.'2) along a position angle of 
146° ± 8°. The size of the continuum source is much larger than the separation of the two 
main protostellar components ( ^0.3"), implying that at least a p art of the emission arises 



from the surrounding envelope. 



Looney. Mundy fc WelchI (119971 ) found that, at 3 mm 



about half of the continuum emission arises from the circumstellar dust disks of the two 
main protostellar components (their observation did not have sufficient angular resolution 
to spatially separate the candidate third protostellar component), and about half from the 
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surrounding dust envelope. The continuum source detected here has an orientati o n sim ilar 
to but size much smaller than the pseudodisk traced in C^^O by 



Momose et al 



([19981) 



suggesting that the extended part of the continuum emission traces the inner regions of the 
pseudodisk. 

Both free-free emission from ionized jets and thermal emission from dust can contribute 
to the continuum emission. In L1551 IRS5, as in many protostellar systems, the ionized jets 
dominate the emission at cm wavelength, whereas dust contributes an increasingly larger 
fraction of the emission towards shorter wavelengths. To assess the relative contribution 
from each component at a given wavelength, we have collected all relevant measurements 
of the continuum emission for L1551 IRS5 at centimeter to submillimeter wavelength as 
listed in Table 2. Figure 2 shows the spectral energy distribution (SED) of the continuum 
emission based on the tabulated data. We modeled the SED as the sum of two power-law 
components with each component having F^^ oc z/", where is the intensity at a given 
frequency v and a the power-law index. A chi-squared fit to the data gives one component 
with a = —0.1 ±0.1 that dominates at low frequencies (short-dashed line in Fig. 2), 
consistent with optically-thin free-free emission. The other component has a = 3.0 ±0.1 
that dominates at high frequencies (long-dashed line in Fig. 2), consistent with thermal 
emission from dust. As can be seen, the contribution from free- free emission to the 
continuum at 1.3 mm is only ~0.1%. 

Assuming that the dust temperature (T^) is uniform, we can estimate the total dust 
mass of the continuum source from the relationship 

M = — p^4^ , (1) 
Bu{Td) K,u 

where d is the distance to the source, B^, is the blackbody function, the dust temperature 
and the dust mass opa city coefficient given by Kj, = 0.1 (z//10^^ Hz)'^ cm^ g~^ 



(IBeckwith &: Sargent 



199ll ). The quantity P is related to the spectral index, a, of 
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the dust continuum emission by the relationship a = 2 + (3. Assuming = 47 K 



( iMoriarty-Schieven et al 



19941 ). the dust mass inferred from Equation (1) is 0.06 ±0.01 



©■ 



This value is in r easonable agreement with those derived at 2.7 mm that are sensitive to the 



same size scales (lOhashi et al. 



1991 



Momose et al. 



19981 ). 



3.2. CO (2-1) 

Figure 3 shows our CO (2-1) channel maps of L1551 IRS5. For presentation purposes, 
none of the maps displayed in this manuscript have been corrected for the primary 
beam response of the SMA; the quantities we computed were, of course, derived from 
maps corrected for the primary beam response. Emission is detectable over the velocity 
range 0.2-10.5 km s~^ measured with respect to the local standard of rest (Vlsr). We 
estimate based on the symmetry of the features observed (see below) a systemic velocity 
of ~6.5 km s~^ for th e L1551 IRS5 protostellar system, comparable to previous estimates 
(IFridlund et al.ll2002l ). The absence of emission within ~0.5 km s"^ of the systemic velocity 
is likely caused by extended emission (mostly from the ambient molecular cloud) that is 
resolved out by the interferometer. 



3.2.1. X-shaped structure 

Figure 4 shows the integrated CO (2-1) intensity map at blueshifted velocities spanning 
0.2-5.5 km s~^ (blue contours) and redshifted velocities spanning 7.4-10.5 km s~^ (red 
contours) separately. To make this map, we first applied Hanning smoothing (a triangular 
function of width five channels at zero intensity) to the channel maps, and then excluded all 
features in the channel maps where the intensity in the corresponding smoothed maps is less 
than 1.5cT. The most prominent feature in this map is an X-shaped structure (as indicated 
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by the diagonal lines) that has its blueshifted pair of arms lying west and its redshifted 
pair of arms lying east of the protostellar system. In the channel maps of Figure 3, the 
blueshifted part of this X-shaped structure can be most clearly seen at velocities spanning 
2.3-5.2 km s~^ and the redshifted part 7.6-10.5 km s~^. Both pairs of arms have an opening 
angle of ~90°, and extend ~20"-30" (^3000-4000 AU) from center. They have a symmetry 
axis along the NE-SW direction at a position angle of about 70°, which is close to the 
orientation of the major axes of the ionized jets as well as the large-scale bipolar molecular 
outflow. Furthermore, the X-shaped structure exhibits the same velocity pattern (i.e., 
redshifted to the E and blueshifted to the W) as the large-scale bipolar molecular outflow. 

In both the channel maps of Figure 3 and integrated-intensity map of Figure 4, the 
blueshifted NW arm can be seen to be much broader than the other three arms of the 
X-shaped structure. A close examination of the channel maps reveals that the NW arm 
comprises several distinct and relatively compact (bullet-like) spatial-kinematic components. 
For example, the broadest portion of the NW arm at about 10" from center (see Fig. 4) 
is dominated by a pair of comparatively high-velocity bullets that can be best seen in the 
channel maps (Fig. 3) at velocities of 2.1 km s~^-3.4 km s~^. The tip of the NW arm 
is dominated by two lower-velocity bullets that can be best seen in the channel maps at 
3.7 km s~^-4.7 km s^^. On the other hand, the two redshifted arms each comprise only 
one spatial- kinematic component (see below). In this way, the properties of the X-shaped 
component seem to mirror that of the large-scale bipolar molecular outflow. Several distinct 
spatial-kinematic components have been seen in single-dish CO (2-1) obser vations spaced 



at se mi-regular intervals along the blueshifted outflow lobe of L1551 IRS5 (Bachiller et al. 



19941 ). Curiously, no such distinct components are present in the same observations of the 
redshifted outflow lobe. 

Figure 5 shows position-velocity (PV-) diagrams along each opposing set of arms as 
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indicated by the diagonal lines in Figure 4. The straight hnes drawn in Figure 5 are not 
meant to guide the eye, but instead are our model fits to the PV-diagrams as we shall 
describe later in §4.1. The PV-diagrams along the NE (Fig. 5(a), upper right quadrant) and 
SE (Fig. 5(b), upper right quadrant) arms look very similar, with each having a velocity at 
peak intensity that increases with radial distance (i.e., Hubble-like expansion) out to ~20" 
from center. The SW arm (Fig. 5(a), lower left quadrant) is only weakly detected in the 
PV-diagram, but nevertheless can be seen to have a velocity gradient comparable to the NE 
and SE arms. On the other hand, the PV-diagram along the NW arm (Fig. 5(b), lower left 
quadrant), shows at least two distinct velocity components. The component labeled NW-a 
increases comparatively quickly in velocity with radial distance out to ~10" from center. 
The component labeled NW-b increases more slowly in velocity with radial distance out to 
~30" from center. NW-a corresponds to the pair of higher-velocity bullets, and NW-b to 
the pair of lower-velocity bullets, mentioned above. Neither of these components have a 
velocity gradient that matches the comparable velocity gradients of the other three arms of 
the X-shaped component. 

3.2.2. S-shaped structure 

The next most prominent feature in the integrated-intensity map of Figure 4 is located 
between the blueshifted and rcdshifted arms of the X-shaped structure and extends to 
~10" from center. The outer portion of this feature on the NE side appears to twist north 
and that on the SW side twist south, creating a S-shaped structure. Its symmetry axis is 
oriented at a position angle of about 55°, approximately aligned with the symmetry axis of 
the X-shaped component as well as major axis of the large-scale bipolar molecular outflow. 

The S-shaped structure exhibits blueshifted velocities NE of center, and redshifted 
velocities SW of center. The velocity pattern of the S-shaped structure is therefore opposite 
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to that of both the X-shaped structure and large-scale bipolar molecular outflow. In the 
channel maps of Figure 3, the blueshifted part of the S-shaped structure can be seen most 
clearly at velocities spanning 3.4-5.5 km s~^ and the redshifted part 7.4-8.2 km s^^. The 
line-of-sight velocity of the S-shaped structure therefore does not reach velocities as high 
(far from Vsys) as the X-shaped structure. As described in §4.3.21 this structure exhibits an 
approximately sinusoidal variation in its radial velocity along its length. 

3.2.3. Compact central component 

Finally, a compact component centered on the location of L1551 IRS5 can be seen in 
the channel maps of Figure 3 extending to higher blueshifted and redshifted velocities than 
either the X- or S-shaped components. It is the only detectable or dominant component 
visible at blueshifted velocities spanning 0.2-3.4 km s^^ and redshifted velocities spanning 
9.2-10.5 km s~^. In the integrated-intensity map of Figure 4, although difficult to visually 
separate from the base of the X- or S-shaped structures, the compact central component 
comprises the closest redshifted and blueshifted intensity peaks straddling the L1551 IRS5 
protostellar system. 

To better extract the properties of the compact central component, we have fitted a 
two-dimensional gaussian structure to this component at each velocity where it is detectable 
in the channel maps. Figure 6 shows the integrated C0(2-l) intensity map made from the 
derived fits with the blueshifted and redshifted emission plotted separately. The centroids 
of the redshifted and blueshifted emissions are displaced by ~1.4" (~200 AU) NE and 
SW from center, respectively, at a position angle of about 65°. The kinematic axis of this 
component is therefore closely aligned with the symmetry/major axis of, and has the same 
velocity pattern as, the X-shaped component and large-scale bipolar molecular outflows. 
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3.2.4- Base of X- shaped and S-shaped components 

To better study the structure of the X- and S-shaped components near their base, 
we have subtracted the compact central component from the channel maps using the 
above-mentioned fits to this component ( §3.2.3I) . Figure 7 shows the integrated C0(2-l) 
intensity maps made before (upper row) and after (lower row) the compact central 
component was subtracted, separated into two different velocity regimes. The higher 
velocity regime (Fig. 7, left column) spans 0.2-3.4 km s~^ (i.e., < Vlsr — 3.1 km s~^) 
at blueshifted velocities (blue contours) and 8.9-10.5 km s~^ (> l^sR + 2.4 km s~^) at 
redshifted velocities (red contours). The lower velocity regime (Fig. 7, right column) 
spans 3.7-5.5 km s"^ (> Vlsr — 2.8 km s~^) at blueshifted velocities (blue contours) and 
7.4-8.7 km s^^ (< Vlsr + 2.2 km s~^) at redshifted velocities (red contours). Like before, 
to make these maps we first applied Hanning smoothing (a triangular function of width 
five channels at zero intensity) to the channel maps, and then excluded all features in the 
channel maps where the intensity in the corresponding smoothed maps is less than 1.5cr. 

The S-shaped component, which as mentioned above does not reach velocities as 
high as the X-shaped component, is barely detectable in the higher velocity regime. In 
this regime, the compact central component can be well separated from the X-shaped 
component as can be seen by comparing the map made before and after subtraction. In the 
lower velocity regime, all three components are detectable. The compact central component 
dominates the emission close to the protostellar system as can be seen by comparing the 
map made before and after subtraction. 
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3.2.5. Physical parameters 

We have estimated physical parameters for all three components separately after 
correcting our maps for the primary beam response of the SMA antennas. For simplicity, we 
assumed that the observed CO (2-1) emission is optically thin, in local thermal equi l ibrium 



(LTE), and has an excitation temperature of 10 K as derived by IStojimirovic et al.l (120061 ) 
based on the intensity ratio between CO (3-2) and CO (1-0). We also assumed a H2-to-CO 
abundance ratio of 10^. We adopted an inclination of 30° with respect to the plane 
of the sky for all three components, such that their major axes are orthogonal to the 
equatorial plane of the circumstellar disks of the two main protos tellar components and 



their surrounding flattened condensation ( ILim fc Takakuwa 



20061 ). As we shall discuss 



later (§4.3), the adopted inclination is probably appropriate for the X-shaped and compact 
central components, and close to median value in the allowed range from our model fit for 
the S-shaped component (but with the NE side the near side, opposite to the X-shaped and 
compact central components). 

In Table 3, we tabulate for each component its estimated mass M, outflow momentum 
P = J2v ^(^) kinetic energy E = \ Ylv ^iy) mass loss rate M = J2v 
dynamical age corresponding to the maximum derived dynamical time t{v) = v/r, where r 
is the distance from the driving source to the emission peak of the outflowing gas at a given 
velocity v. Note that all these values except for dynamical age are lower limits, in part 
because we have assumed that the CO (2-1) emission is optically thin, and in part because 
emission close to the systemic velocity was entirely resolved out. The X-shaped component, 
which is the most obvious component in our map, dominates in mass, momentum, and 
energy. It also has the oldest dynamical age. The mass-loss rate is highest for the compact 
central component, which also has the youngest dynamical age. 



4. Discussion 



Our C0(2-l) map at an angular resolution of 2'/3 x 4'.'7 reveals, for the first time, 
three distinct bipolar molecular outflow components in the vicinity (< 4000 AU) of the 
L1551 IRS5 protostellar system. In the following, we shall discuss the likely nature of each 
of these components as well as their driving sources. 



4.1. X-Shaped Component 



The X-shaped component has a symmetry axis closely aligned with and exhibits the 
same velocity pattern as the large-scale bipolar molecular outflow. It has the largest mass, 
momentum, and energy of the three outflow components that we detected in the vicinity of 
the L1551 IRS5 protostellar system. This X-shaped component resemb les similar structures 



seen around a number of other pro tostars such as IRSl i n Barnard 5 (IVelusamy fc Langer 



998h. IRAS 20582+77 24 in L1228 flArce fc Sargent 



(lArce &: Sargent 



20041 ) ■ and IRAS 5295+1247 in L1528B 



20051 ). In all these cases, the X-shaped structure observed is thought to 



comprise the limb-brightened walls of a cone-shaped cavity entrained and excavated by 
outflows from the central protostars. Indeed, scattere d light from a cone-shaped cavity can 



be seen in the near-IR on the SW side of L1551 IRS5 (jHodapp 



19941 ) with an opening angle 



comparable to that measured for the (blueshifted) arms of the X-shaped component. (The 
NE side of the cavity is presumably obscured in the near-IR by material in the condensation 
and extended envelope.) 



Arce fc Sargent 



(120061 ) have conducted a ^^CO survey of molecular outflows within 
<10'^ AU of nine protostellar (comprising a mixture of both single and multiple) systems 
at different evolutionary stages. They found a clear difference in the morphology of the 
molecular outflow lobes between Class and more evolved Class I protostars. Class 
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protostars exhibit relatively collimated jet-like outflows or cone-shaped lobes with opening 
angles < 55°, whereas Class 1 protostars have outflow lobes with wider opening angles of 
> 75° (with ~125° the largest openin g angle seen). L1551 IR S5 is classifled (based on its 



Infrared SED) as a Class I protostar (IKeene fc Masson 



1990l ). and so the opening angle of 



its outflow cavity (~90°) is comparable with that seen around other protostars at a similar 
evolutionary stage. 



Pyo et al. 



(12002 



20051 ) have revealed two 



Near-IR observations in the [Fe II] line by 
velocity components in the ionized outflow from the main northern protostellar component 
in L1551 IRS5. One of these outflow components comprises a high-velocity (~300 km s~^) 
and narrowly-coUimated wind, and the other a low-velocity (~100 km s~^) and wide-angle 
wind. The latter has an opening angle of ~100°, which is similar to the opening angle of 
the cavity seen through scattered light in the near-IR and as the X-shaped component. The 
X-shaped component may therefore have formed in two steps: flrst the high-velocity and 



narrow. 



Stabler 



^-col 



imated wind passes through and in the process entrains molecular gas (e.g.. 



19941 ). followed by a low- velocity and wide-angle wind that pushes for ward and 



Bachiller et al. 



(119941 ) have found 



aside the molecular gas to excavate a cone-shaped cavity, 
that the large-scale bipolar molecular outflow associated with L1551 IRS5 exhibits several 
successive high-velocity features along its blueshifted lobe suggestive of different ejection 
episodes. The X-shaped component may therefore comprise (primarily) the base of that 
portion of the larger-scale outflow created in the latest ejection episode. 

We now investigate what properties a cone-shaped outflow cavity should have to 
reproduce both the observed structure and kinematics of the X-shaped component. Figure 8 
shows a sketch of a cone with a half-opening angle 9 and its symmetry axis inclined by 
an angle i to the plane of the sky. As suggested by the PV-diagrams of Figure 5, we 
assume that the cone exhibits a Hubble-like radial expansion such that its radial velocity 
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is given by Vr = Cr, where C is the coefficient of the radial flow and r the radial distance 
from the driving source. We assume i = 30° so that the symmetry a xis is orthogonal t o 



the c ircumstellar disks of the two main protostellar components (see 



Lim fc Takakuwa 



20061 ). and measure from the integrated- intensity map of Figure 4 a half-opening angle of 
6 ^ 45°. The only free parameter is therefore C, which in the PV-diagram of Figure 5 
specifies the observed velocity gradient along the arms of the X-shaped component. We 
found that a satisfactory fit to all except the blueshifted NW arm can be obtained with 
C = 0.18 km s~^ arcsec"^, and plot the model fit as solid lines in the PV-diagram of 
Figure 5. As mentioned earlier in §3.2.1, the blueshifted NW arm exhibits two distinct 
velocity components in the PV-diagram of Figure 5, neither of which have velocity gradients 
that matches the otherwise similar velocity gradients along the other three arms. 

The acceleration constant of the Hubble-like radial expansion estimated for the 
X-shaped component of C ~ 0.18 km s~^ arcsec"^ ~ 250 km s~^ pc~^ is much higher than 
that estimated from single-dish observ ations for the large-s cale bipolar molecular outfiow 
from L1551 IRS5 of ~70 km s^^ pc^^ (iFridlund et al.lll989h . Interferometric observations 



of the bipolar molecul ar outfiow 



of ~490 km s ^ pc ^ (lYen et al. 



Tom B 335 reveal an even higher acceleration constant 



20091 ). On the other hand, the acceleration constant 



of the Hubble-like radial expansion of large-scale (> 0.1 pc) molecular o utfiows in other 



low-m ass protostars have measured val ues of ~65 km s ^ pc ^ in L1157 (IBachiller et al. 



2001), ~97 km s"^ pc"^ in both L483 flTafalla et al. 



19971 ). and ~120 km s"^ pc"^ in IRAM 04191+1522 (lAndre et al. 



2000h and BHR 7 1 flBourke et al 



19991 ). All these values 



are lower than that measured in interferometric observations of either L1551 IRS5 or B335 
that probe the base of their respective outfiows, suggesting that the acceleration constant 
of the Hubble-like radial expansion decreases outwards along a given outfiow lobe. 
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4.2. Compact Central Component 



The compact central component has its major axis aligned with the symmetry axis 
of the X-shaped component, as well as major axis of the large-scale bipolar molecular 
outflow. Furthermore, it exhibits the same velocity pattern as both these more spatially 
extended components. Together with the fact that it reaches much higher blueshifted and 
redshifted velocities than the X-shaped component, the compact central component most 
likely comprises material lying along the symmetry axis of the cone-shaped cavity that 
has been entrained relatively recently by the jet(s) from one (or both) of the two main 
protostellar components. High angular-resolution observations of the bipolar molecular 
outflows from other protostars also sometimes show a relatively collim ated outflow lying 



along the symmetry axis of a cone-shaped outflow cavity (e.g., HH211; 



Gueth fc Guilloteau 



19991 ). with the more collimated component attributed to molecular gas entrained by a jet. 



Fridlund fc Liseaul (119981 ) found that the two optical jets visi ble on the S W side of 



LI 551 IRS5 have velocities of ~300 km s ^. As mentioned in 54.1, 



Pyo et al. 



(12002 



20051) 



found that the comp onent at higher velocities of ~300 km s ^ in the northern optical jet is 



narrowly collimated. iFridlund fc Liseaul (119981 ) inferred based on the observed Ha and [S 11] 
line intensities a particle density of ^10"^ cm~^ in these jets, and therefore a momentum 
for each jet of ~7 x 10~^ km s~^. The combined momentum of the optical jets is 
therefore comparable with the momentum we inferred for the compact central component 
of ~1.33 X 10^^ Mq km s^^, demonstrating that the optical jets indeed have sufficient 
momentum to drive the compact central component. At their velocities of ~300 km s~^, 
the optical jets would travel the entire length of the compact central component in just 
~30 yrs. The presence of newly entrained molecular gas in the outflow cavity may indicate 
that a significant quantity of molecular gas has recently been injected into this cavity from 
the surrounding infalling envelope, which tries to fill or is diverted into the cavity at or near 
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its base. Alternatively, or in addition, a relatively recent enhancement in the activity (i.e., 
mass-loss rates) of one or both optical jets have led to an enhancement in the amount of 
entrained material in the cavity. 



Using the Chandra X-ray Observatory, iBally et al.l (120031 ) have detected a compact, 
but slightly resolved. X-ray source located about 0.5"-l" south-west of L1551 IRS5 at the 
base of the two optical jets. They proposed several models for how the observed X-ray 
emission may be produced, all of which involve radiating plasma downstream of a shock 
where either: i) a wide-angle flow, launched and accelerated within a few AU of one member 
of the IRS 5 protobinary, is redirected and collimated into a jet; ii) one or both jets collide 
with a circumbinary disk; and (iii) the two jets collide. Here, we propose that shocks 
produced when one or both jets slam into molecular gas diverted into the X-shaped cavity 
provide a simpler if not more plausible explanation for the observed X-ray emission. The 
resulting entrained molecular gas would be expected to extend beyond the location of the 
X-ray source, as is observed for the compac t central conapone nt. To produce the observed 



X-ray luminosity and plasma temperature. 



Bally et al. 



(120031 ) find that a preshock gas 



density of 1-10 x 10^ cm~'^ and shock velocities larger than 350 km s^^ are required. Such 
molecular hydrogen gas densities are required to excite the observed CO (2-1) emission 
from the compact central concentration. The line-of-sight velocity of at least one of the two 
op tical iets reach ^3 km s'^, and if inclined by ~30° to the plane of the sky as inferred 



by 



Lim fc Takakuwal (120061 ) would therefore reach true velocities of ~600 km s ^. Such 



high jet velocities can produce the required shock velocities. 

Instead of tracing an outflow, could the compact central component constitute the 
inner region of the flattened molecular condensation around L1551 IRS5? In such a case, 
the observed NE-SW velocity grad ient would then correspond to infalling motion as seen 



in C^^O(l-O) (iMomose et al 



19981 ). The measured velocity gradient caused by the infalling 
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motion at the same scale as the com pact CO (2-1) centra. 



■^0.13 km s ^ arcsec ^ (see Fig. 3 in 



Momose et al 



component (~300 AU) is at most 



19981 ). For comparison, the measured 



velocity gradient of the co mpact central compone nt is at least ~0.60 km s ^ arcsec ^ 



Furthermore, as shown by 



Takakuwa et al. 



(120041 ). the inner region of the condensation 



is dominated by rotation and not infall, and exhibits a velocity gradient in the NW-SE 
direction (along the major axis of the condensation) and not NE-SW direction. Thus, the 
compact central component is unlikely to comprise the inner regions of the condensation 
around L1551 IRS5. 



4.3. S-shaped Component 

The S-shaped component would not have been spatially resolved and hence recognized 
in previous single-dish CO observations. Given that it exhibits a velocity pattern opposite 
to the compact central component, X-shaped component, and large-scale bipolar molecular 
outflow, at first sight it is unlikely to comprise material entrained by the collimated winds 
of the two main protostellar components. Furthermore, this component exhibits the lowest 
(projected) maximum detectable velocity of all three components. Below, we consider three 
possible explanations for this component: i) the near and far sides of the outflow cavity 
walls; ii) a precessing outflow; and iii) an outflow that winds in space due to the orbital 
motion of its driving source. 



4-.3.1. Near and Far Sides of Outflow Cavity Walls 

Here, we consider the possibility that the S-shaped component is part of the same 
outflow cavity traced by the X-shaped component. In the model described in §4.1, the 
outflow cavity has a half-opening angle (~45°) that is larger than its inclination to the plane 
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of the sky (~30°). On the NE side of the outflow cavity where the observed hmb-brightened 
waUs are redshifted, the waUs on the near side of this cavity is therefore blueshifted. 
Conversely, on the SW side of the outflow cavity where the observed hmb-brightened 
waUs are blueshifted, the walls on the far side of this cavity is redshifted. The S-shaped 
component may therefore comprise the near side of the outflow cavity on the NE side and 
the far side of the outflow cavity on the SW side. 

In such a case, the model described in §4.1 that flts the observed morphology 
and kinematics of the X-shaped component should also fit the observed morphology 
and kinematics of the S-shaped component. Figure 9 shows the PV-diagram along the 
symmetry axis (position angle of 70°) of the X-shaped component. Emission close to 
the center straddled by the horizontal short-dashed lines is dominated by the compact 
central component. The solid line shows what our model described in §4.1 (with 
C = 0.18 km s~^ arcsec"^ and i — 30°) predicts if the S-shaped component comprises the 
near side of the outflow cavity on the NE side and the far side of the outflow cavity on 
the SW side as described above. Clearly the S-shaped component cannot be fit by such 
a model. The PV-diagram of the S-shaped component in Figure 9 is better fit by the 
long-dashed line, which corresponds to an acceleration constant for a Hubble-like radial 
expansion of C = 0.75 km s~^ arcsec"^ for the same assumed inchnation of i = 30°. The 
long- dashed line in Figure 5 shows what this model predicts for the PV-diagram of the 
X-shaped component. Clearly, the X-shaped component cannot be flt by such a model. 

To simultaneously fit the PV-diagrams of both the X- and S-shaped components in 
a single model for the cone-shaped cavity, this cavity is required to have an inclination 
angle of i = 16° and acceleration constant for its Hubble-like radial expansion of 
C = 0.36 km s^^ arcsec^^. With this set of parameters, the model fit to the X-shaped 
component is closely indicated by the solid line in Figure 5 and to the S-shaped component 
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the dashed line in Figure 9. The symmetry axis of such a cone-shaped cavity would then be 
tilted by ~14° from a direction orthog onal to the circumstellar disks of the main northern 



and south ern protostella.r cona ponents (ILim &: Takakuwa 



envelope (IMomose et al. 



20061 ) as well as the surrounding 



19981 ). This model does not provide an obvious explanation for 
why the far wall on the NE side and near wall on the SW side are not detected. More 
problematically, it does not provide a natural explanation for the S-shaped morphology of 



this component nor its oscillatory line-of-sight velocity as described next. 



4-3.2. A Precessing Outflow 

The morphology of the S-shaped component resembles a precessing outflow. To assess 
whether its kinematics also resemble a precessing outflow, we adopt the following simple 
model. We assume that the outflow is driven by a jet with a constant velocity Vj^t that 
precesses at a constant angular frequency with a period Tp. Such a jet will describe a 
sinusoidal pattern in the sky with a deflectio n amplitude that inc reases with distance from 



the driving source in the form (adapted from lEisloffel et al.lll996l ): 

cosip —sinip \ I a • / • sin(27r//A + 0o) 



(2) 

sinip cosip j \ I ■ cos i 

where x and y are the usual Cartesian coordinates, a the precession amphtude, A the 
precession length scale, / the distance from the source, 0o the initial phase at the source, ip 
the position angle of the outflow symmetry (precession) axis in the plane of the sky, and i 
the inclination angle of the outflow symmetry axis to the plane of the sky. The line-of-sight 
velocity, li^os, can be derived from the equation: 

^L05 = ^LSfi ± Vc-o [cos 6^ sin z -|- sin 6' cos 2 cos(27r//A + 0o)] , (3) 

where Vco is the measured C0(2-l) velocity, and 6 the half-opening angle of the winding 
outflow. Figure 10 shows a schematic picture of this model in which the jet precesses 
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clockwise about its symmetry (precession) axis. Note that smaller inclinations (i) require 
smaller precession amplitudes {a) to produce a given line-of-sight velocity (Vlos) deflection 
in the sky. 

We start by considering a straight jet. To fit the observed morphology and kinematics, 
we find that the inclination of the jet symmetry axis to the plane of the sky is restricted 
to the range 15° < i < 45°. For illustration, we take i ~ 30°, allowing us to directly 
compare the energetics of the S-shaped component with the X-shaped and compact central 
components (see §3.2.5). For this inclination, we find that the model parameters listed in 
Table 4 (first row) provide the closest fit to both the observed morphology and kinematics 
of the S-shaped component. Within the allowed inclination range, other sets of model 
parameters can provide a more or less equally satisfactory fit to the observed morphology 
and kinematics. 

Figure 11 compares the model predictions (black line) in the case where i ~ 30° 
with the observed morphology (upper panels) and line-of-sight velocity (lower panel) of 
the S-shaped component. Here, the jet that drives the S-shaped component precesses at 
an angle oi 9 ^ 22° to its symmetry (precession) axis. The morphology of the S-shaped 
component is well reproduced in this model where the driving jet lies outside the conical 
outflow cavity traced by the X-shaped component along most of its length, thus enabling it 
to entrain the surrounding material. The line-of-sight velocities plotted in Figure 11 (lower 
panel) are derived by fitting the measured spectra at regularly spaced intervals along the 
predicted path of the precessing outfiow with simple gaussians. Two sets of error bars 
are plotted for each point, the smaller (solid line) corresponding to the uncertainty in the 
measurement of the central velocity, and the larger (dashed fine) the FWHM of the fitted 
gaussian. As can be seen, the line-of-sight velocity also changes more or less sinusoidally in 
reasonable — albeit not perfect — agreement with the model. Given that the entrained 
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material may have had different initial velocities along the length of the outflow reflecting 
the rotation and contraction of the surrounding envelope, significant differences between 
the predicted and measured velocities may not be unexpected. 

The precession period, Tp, can be determined from the wavelength A and jet velocity 
Vjet such that 



(y.s) = 4.76 X , (4) 



where d is the distance to L1551 IRS5 (140 pc). Assuming a typical jet velocity of 
~50-150 km s~^ (see below for reason why these values were adopted), the precession 
period is ~130-400 yrs. 

To see if we can produce an even better fit to the observed morphology and, especially, 
kinematics, of the S-shaped component, we relax the requirement that the jet driving this 
component be straight. A bent jet (see inset in Fig. 10) — in which the two opposing arms 
are not aligned — have been seen or invoked to explain the bipolar molecular outflows from 
a number of prot ostellar systems. The me chanisms that can lead to such bending have 



been explored by iFendt fc Zinneckerl (Il998l ): in the present situation, wide-angle winds at 
the base of the main northern and/or southern protostellar components pushing on the 
jet of the third protostellar component provides yet another possibility. Adopting once 
again i ~ 30°, we find that the model parameters listed in Table 4 (second row) provide a 
better fit to the observed kinematics. Figure 12 plots the model predictions (black line) 
for the observed morphology (upper panels) and line-of-sight velocity (lower panel) of the 
S-shaped component, where the jet that drives this component precesses at an angle of 
~10° to its symmetry axis. As can be seen, the measured line-of-sight velocities along the 
precessing outflow are in relatively good agreement with model predictions for a bent jet 
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where one arm lies at an angle of 163° to the other (for a straight jet, the two arms lie at an 
angle of 180°). For a jet velocity of ~100 km s~^ as assumed above, the precession period 
corresponds to ~135 yrs. 

Could the S-shaped component, which as we have shown can be well fitted by a 
precessing out flow, be driven by the ca ndidate third protostellar component in L1551 IRS5 



discovered by iLim fc Takakuwal (l2006l )? The circumstellar disk of this third component is 
tilted with respect to the closely aligned circumstellar disks of the two main protostellar 
components, as well as their surrounding pseudodisk. The outflow driven by a jet from 
the third protostellar component should therefore have a different inclination and/or 
position angle in the sky compared with the outflows driven by the northern and southern 
protostellar components. Indeed, the S-shaped outflow is tilted in the opposite sense of the 
sky with respect to the compact central outflow, X-shaped outflow cavity, and large-scale 
bipolar molecular outflow, all of which are driven by the two main protostellar components. 
In the model incorporating a straight jet, the jet currently emerges at its origin at a position 
angle of 78° in the plane of the sky. By comparison, the major axis of the circumstellar 
disk of the third component is at a position angle of 118° ± 8°, and hence the jet is bent by 
~50° from orthogonal to the disk. In the model incorporating a bent jet, at or close to its 
origin the eastern (blueshifted) arm of the jet emerges at a position angle of ~65° and the 
western (redshifted) arm of the jet at a position angle of ~48° as sketched in the insert of 
Figure 10, within ~20° and ~37° respectively from orthogonal to the disk. 

Tidal interactions between the third protostellar component and its more massive 
neighbor(s) c ould cause t he dis k of this component to precess, resulting in the precessing 



outflow seen. 



Bate et al. 



(120001 ) have shown that such interactions between two protostellar 
components in a binary system should cause their circumstellar disks to precess at a 
frequency that is a small fraction (of order a few percent) of their orbital angular frequency. 
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Taking a mass for the main northern protostellar component of 0.45 Mq and the third 
protostellar component of 0.1 Mq, and a physical se paration between the tw o components 



Lim fc Takakuwa 



20061 ). the orbital 



equal to their projected separation of ~13 AU (see 
period of the third protostellar component about the northern protostellar component is 
~63 yrs. For comparison, the estimated precession period of the S-shaped outflow is a 
few hundred years, assuming a jet velocity for the third protostellar component that is 
proportional to the square- root of its mass (i.e., proportional to the Keplerian velocity). 
The results are compatible with the idea that the circumstellar disk of the third protostellar 
component is precessing at an angular frequency much smaller than the orbital angular 



frequency of this component about the main northern protostellar component. 



4 .3. 3. A Winding Outflow due to Orbital Motion 



As mentioned by iFendt &: Zinneckerl (119981 ) . the orbital motion of a protostar about 
its binary companion could also cause its originally straight jet to be deflected backwards. 
Because such a bent jet possesses mirror rather than point symmetry (as in the case 



of a precessing iet 



Masciadri fc Raga 



, it s hould drive a winding outflow with reflection symmetry (e.g.. 



2OO2I ). We have therefore examined whether the morphology and 



kinematics of the S-shaped component can be flt by a winding mirror-symmetric outflow. 

For simplicity, we consider a circular orbit. To produce the observed morphology, the 
orbital velocity, Vorb, of the protostar that drives the jet is required to satisfy the relation: 



Vorb COS J 



Vjet COS I 



tan^ 



(5) 



where as before Vjet is the jet velocity, i the inclination to the plane of the sky at which the 
jet emerges from its driving source (i.e., tangent to the jet at its origin), 6 the half-opening 
angle of the winding outflow (also equal to the angle at which the jet is bent from its 
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original trajectory), and j the orbital inclination to the plane of the sky. The resulting 
hne-of-sight velocity of such a jet will oscillate sinusoidally about the jet velocity, with the 
oscillation caused by the orbital motion of the protostar driving the jet. As can be seen 
from Figure 12, the molecular outflow oscillates with a half amplitude of ~0.5 km about 
a line-of-sight velocity of ~1.5 km s~^; i.e., ratio between these two values of ~0.3. To 
produce the observed kinematics, the orbital velocity of the protostar is therefore required 
to also satisfy the relation: 

'-^^^ ^ 0.3 . (6) 

Vjet sm i 

We assume throughout that the orbital diameter of the protostar is much smaller than the 
observable width of the winding outflow (i.e., that the jet emerges from a stationary point), 
an assumption that we shall return to later. 

We found that a mirror-symmetric outflow model can provide a reasonable fit to the 
observed morphology and kinematics of the S-shaped component over only a relatively 
narrow range in values of 14° < ^ < 19° and 40° < i < 50°. For illustration, we take median 
values of ^ ~ 17° and i ~ 45°, in which case Equations (5) and (6) can be solved to give 
j Ri 45°. To fit the observed morphology and kinematics, we also require A = 17". The 
model fit for these values is shown in Figure 13. It is superior, especially in reproducing the 
observed kinematics, to the model fit for a straight precessing jet, but inferior to that for a 
bent precessing jet. 

The orbital period of the protostar is equal to the period of the winding outflow, 
which has the same dependence on the model parameters as the precession period given 
by Equation (4). Substituting into this equation Vjet ~ i'or6/0.3 from Equation (6) given 
the model parameters used here, we find Torb{yrs) = 0.3 x 4.76 x A x d{pc)/vorb- The 
product Vorb X Torb is simply the circumference of the orbit, from which we derive an orbital 
diameter of ~236 AU (~1'.'7). Thus, the orbital diameter is indeed much smaller than the 
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observable width of the winding outflow, satisfying an underlying assumption of our model. 
In this model, the minimum projected orbital radius is ~83 AU (~0'.'6), many times larger 
than the observed projected separation of the 3rd component from the main northern 
protostellar component of ~13 AU (~0'.'l). Indeed, the minimum projected orbital radius 
is a factor of ~2 larger than the observed projected separation of the 3rd component from 
the main southern protostellar component. Thus, to fit all the observational constraints, 
elliptical orbits with high eccentricities would have to be considered (or a fourth yet unseen 
protostellar component invoked). 

For the model parameters used here, the orbital period can also be expressed as 

1290(yrs) 



where Mat and Msrd are the masses of the main northern and 3rd protostel 



ar co mponents 



respectively. Taking Mn ^ 0.45 Mq and M^rd ~ 0.1 Mq (iLim fc Takakuwall2006l ). we find 
Torb ~ 1700 yrs and hence Vorb ~ 2.1 km s~^ and Vjet ~ 7 km s~^. The jet velocity in this 
model is unreasonably small, a serious flaw in the picture where the winding outflow is due 
to the orbital motion of the driving source. 



5. Summary 

L1551 IRS5 exhibits the first recognized bipolar molecular outflow from a protostellar 
system. Single-dish maps show that this outflow can be traced in both CO (1-0) and 
CO (2-1) to 1.4 pc on the sky along the NE (redshifted lobe) and SW (blueshifted lobe) 
sides of the system along a position angle of ~50°. This outflow is presumably driven by 
the two main protostellar components in L1551 IRS5, each of which exhibits an ionized jet 
that is closely aligned with the major axis of the larger scale molecular outflow. 

Here, we imaged the molecular outflow within ~4000 AU of L1551 IRS in C0(2-l) 
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with the SubMiUimeter Array (SMA) at an angular resolution of ~4" (linear resolution of 
~560 AU). Our observations revealed three distinct spatial-kinematic components: 

1. An X-shaped structure centered on L1551 IRS with arms that can be traced to 
~20"-30" (~2800-4200 AU) from the system. It has a symmetry axis along the NE 
(redshifted side) to SW (blueshifted side) direction at a position angle of ~70°. This 
component is therefore closely aligned with, and exhibits the same velocity pattern 
as, the large-scale bipolar molecular outflow. 

2. A previously unknown S-shaped structure centered on L1551 IRS5 that can be 
traced to ~15" (~2100 AU) from the system. It has a symmetry axis along the 
NE (blueshifted side) to SW (redshifted side) direction at a position angle of ~55°. 
Although also closely aligned with the major axis of the large-scale bipolar molecular 
outflow, this component exhibits a velocity pattern in the opposite sense. There is no 
known counterpart to this outflow on large (> 0.1 pc) scales. This component has the 
smallest linewidth of all three components. 

3. A previously unknown compact component centered on L1551 IRS5 that can be 
traced to ~1.4" (^200 AU) from the system. It has a major axis along the NE 
(redshifted side) to SW (blueshifted side) direction at a position angle of ~65°. This 
component is therefore closely aligned with, and exhibits the same velocity pattern 
as, the large-scale bipolar molecular outflow. It has the largest linewidth of all three 
components. 

We interpreted our results in the following manner: 

1. The X-shaped component comprises the limb-brightened walls of a cone-shaped 
outflow cavity. The molecular gas originally in this cavity was likely entrained by 
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the high-velocity colhmated winds of the two main protostellar components, and 
subsequently excavated by the low-velocity wide-angle winds of the same components. 
The opening angle of this cavity is ~90°, comparable to that seen in other Class I 
protostellar systems. 

2. The S-shaped component is unlikely to be the front and back walls of the same outflow 
cavity whose limb-brightened walls are seen as the X-shaped component. Instead, 
the S-shaped component is most likely a precessing outflow whose symmetry axis is 
inclined in the opposite sense to the plane of the sky compared with the X-shaped 
and compact central components, as well as the large-scale bipolar molecular outflow. 
This precessing outflow m ay driven by the candida te third protostellar component in 



L1551 IRS5 discovered by iLim fc Takakuwal (120061 ) . who found that the circumstellar 



disk of this component is misaligned relative to the two main protostellar components. 
Gravitational interactions between the third protostellar component and its likely 
more massive northern (and perhaps also southern) neighbor (s) may be causing the 
circumstellar disk and hence jet of this component to precess. 

3. The compact central component likely comprises material newly entrained by one or 
both jets from the two main protostellar components. The relatively short dynamical 
age for this component of ~200 yrs may suggest a recent enhancement in the strength 
of one or both ionized jets responsible for driving this outflow. 



Previous single-dish observations sug gest multiple election episodes in the large-scale 



bipolar molecular outflow of L1551 IRS5 (IBachiller et al. 



19941 ). Our study reveals multiple 



bipolar molecular outflow components near the base of this large-scale outflow. It adds 



additional evidence f or a third protoste l 



to that discovered by 



Lim fc Takakuwal (120061). 



ar co mponent in LI 551 IRS5 that may be related 
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Table 1: Observational Parameters 



PARAMETER VALUE 



Right Ascension (J2000) 04^ SI'" 34n4 

Declination (J2000) 18° 08' 05'/l 

Primary beam (FWHM) ~57" 

Synthesized beam (FWHM) 4'/66 x 2'.'30 (P.A. = -66.7°) 

Central frequencies 231.3 GHz, 241.3 GHz 

Frequency (velocity) resolution 203.125 kHz (~0.264 km s~^) 

Gain calibrators 0423-013 (9.4 Jy), 3C 120 (3.3 Jy) 
Absolute flux and passband calibrator Uranus 

rms noise level (continuum) 0.014 Jy beam~^ 

rms noise level (hne) 0.35 Jy beam~^ at 203.125 kHz resolution 
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Table 2: Continuum Flux Density of L1551 IRS5 



\A/ a Arpl PTi erf n 

V V d V ClCli^ uli 


r liiv T)pri tji^AA ( vn \\r\ 


R pfpypTi ppc; 


18 cm 


1.3 ± 0.4 


1 


3.6 cm 


1.48 ± 0.04 


1 


2.0 cm 


2.1 ± 0.2 


1 


1.3 cm 


3.5 ± 0.4 


1 


7 mm 


12.2 ± 1.0 


1 


7 mm 


10.1 ± 0.7 


2 


2.7 mm 


171 ± 19 


3 


2.7 mm 


162 ± 6 


4 


1.3 mm 


1200 


This work 


850 /im 


7230 


4 


450 /xm 


45700 


4 



^ CI) iRodriguez et all (Il998f): ("2) iLim fc Takakuwa 
( 20061); (3) Momose et. al. (1998); (4) 
Moriartv-Schieven et al. ( 20061 ) 
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Table 3: Properties of the Three Outflow Components 







Mass (Af) 
(10-:^ Mq) 


Momentum (P) 
(10-3 M0 km s-i) 


Energy (E) 
(10-*" erg) 


Mass-loss rate (M) 
(10-7 M0 yr-i) 


Dynamical Age 

(yrs) 


X-shaped 


Red lobe 


0.48 


2.28 


12.12 


2.60 


1.3x10"' 




Blue lobe 


0.77 


4.07 


25.25 


3.37 


1.2x10" 


5-shaped 


Red lobe 


0.09 


0.34 


1.45 


0.42 


6.3x103 




Blue lobe 


0.29 


0.94 


3.44 


1.80 


7.4x103 


Central 


Red lobe 


0.29 


0.95 


3.56 


10.98 


2.7x10^ 




Blue lobe 


0.22 


1.33 


9.66 


11.00 


2.0x10^ 
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Table 4: Model parameters for precessing jet 







a 


Q 


A 


00 (red, blue) 


i 


Straight Jet 


57° 


0.40 


22° 


40 


-90°, -90° 


30° 


Bent Jet 


57° 


0.18 


10° 


20 


110°, -135° 


30° 
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34!2 34!0 33!8 33!6 



RA (J2000) 



Fig. 1. — Continuum map of L1551 IRS5 at 1.3 mm. Contour levels are plotted at —2, 2, 10 
20, 30, 40, 50, 60(7, where the rms uncertainty a = 0.014 Jy beam~^. Cr osses indicates the 



position of the two main protostellar components at 7 mm as measured by 



Lim fc Takakuwa 



(120061 ). The synthesized beam is shown as a filled ellipse at the bottom left corner. 



/ 



1 2 3 

log V (GHz) 

Fig. 2. — SED of the continuum emission from L1551 IRS5. The data points and their cor- 
responding ±lcr measurement uncertainties (where large enough to be seen) are tabulated in 
Table 2. The solid curve shows our best chi-squared fit assuming two power-law components 
with each having FyOav"'. One component has a pa —0.1 as indicated by the short-dashed 
line and dominates at low frequencies, and the other a ~ 3.0 as indicated by the long-dashed 
line and dominates at high frequencies (see text) 



-44- 




Fig. 3. — C0(2-l) channel maps of L1551 IRS5. Contours levels are plotted at —3, 3, 6, 9, 
15, 25, 35 and 45(7, where a — 0.35 Jy beam~^. The central Vlsr of each channel is shown at 
the top left corner of each panel. Crosses indicate the positions to the two main protostellar 
components as in Figure 1. The field of view of the SMA at FWHM is shown by the large 
circle, and the synthesized beam by the filled ellipse at the bottom left corner. 
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RA offset (arcscc) 



Fig. 4. — Integrated CO (2-1) intensity map with blue contours showing emission at 
blueshifted velocities spanning Vlsr = 0.23-5.5 km s~^, and red contours emission at red- 
shifted velocities spanning Vlsr = 7.4-10.5 km s^^. For the blueshifted emission, contour 
levels are plotted at —1.3, 1.3, 2.6, 3.9, 6.4, 10.7, 15.0, and 19.3 Jy km s~^ beam~^, while 
for the redshifted emission, contour levels are at —1.0, 1.0, 2.0, 3.0, 5.1, 8.5, 11.9 and 15.3 
Jy km s~^ beam~^. The cross indicates the location of L1551 IRS5 protostellar system. The 
field of view of the SMA at FWHM is shown by the large circle, and the synthesized beam 
by the filled ellipse at the lower left corner. Diagonal lines indicate the orientation of cuts 




2468 10 2468 10 

Fig. 5. — PV-diagrams of the C0(2-l) emission along position angles (a) 25° and (b) 115° 
passing through the arms of the X-shaped component and central protostellar system (see 
Fig. 4). Contour levels are plotted at 0.7, 1.4, 2.1, 2.8, 3.5, 5.5, 7.5, 9.5, and 11.5 Jy beam-^ 
In each panel, the pair of horizontal short-dashed lines straddling the horizontal solid line 
(which indicates the position of the protostellar system) are separated by twice the FWHM 
of the synthesized beam at that particular position angle. Emission within this region is 
dominated by the compact central component (see text). The vertical line indicates the 
systemic velocity that we inferred for L1551 IRS5. The dark solid lines are the predictions of 
our model described in §4.1. The long-dashed hues are the predictions of a model described 
in §4.3.1. Note that panel (b) has been inverted vertically to show that three of the four 
arms (except the blueshifted NW arm) share the same velocity pattern (see text). 
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Fig. 6. — Integrated C0(2-l) intensity map constructed from gaussian fits to the com- 
pact central component in the channel maps (see text). Blue contours shows emission at 
blueshifted velocities spanning 0.23-5.5 km s~^, and red contours emission at redshifted ve- 
locities spanning 7.4-10.5 km s~^. Contours levels are plotted at —3, 3, 9, 15, 25, 35, and 
45(7, where a = 0.43 Jy km s^^ beam^^ for the blueshifted and a = 0.34 Jy km s^^ beam~^ 
for the redshifted emission. The synthesized beam is shown as a filled ellipse at the bottom 
left corner. 




RA offset (arcsec) 



Fig. 7. — Integrated CO (2-1) map with the compact central component included (upper 
row) and with it subtracted (lower row) separated into two velocity regimes (left and right 
columns). In the higher velocity regime (left column), the blue contours span velocities 0.23- 
0.34 km s~^ and the red contours 8.9-11 km s~^. In the lower velocity regime (right column), 
the blue contours span velocities 3.7-5.5 km s~^ and the red contours 7.4-8.7 km s~^. In the 
left column, blue contours levels are plotted at —1.0, 1.0, 2.0, 3.0, 5.1, 8.5, 11.9 and 15.3 
Jy km s^^ beam^^ and red contours levels are plotted at —0.8, 0.8, 1.5, 2.3, 3.8, 6.3, 8.8 and 
11.3 Jy km s~^ beam~^. In the right column, blue contours levels are plotted at —0.8, 0.8, 
1.6, 2.4, 3.9, 6.5, 9.1 and 11.7 Jy km s~^ beam~^ and red contours levels are plotted at —0.7, 
0.7, 1.4, 2.1, 3.5, 5.8, 8.1 and 10.4 Jy km 8"^ heam~\ 
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Sky plane 




Fig. 8. — A schematic picture of a cone-shaped outflow cavity that constitutes our model 
for the X-shaped component. The cone has a half-opening angle of 9, and is inclined by an 
angle i to the plane of the sky. The arrows indicate the direction of gas motion, which are 
along the cavity walls. 
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Fig. 9. — PV-diagram of the C0(2-l) emission along the symmetry axis of the S-shaped 
component along a position angle of = 70°. Contour levels and the different lines are the 
same as in Figure 5. 
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Fig. 10. — A schematic picture of a precessing jet, represented by the blue/red sinusoidal 
curve, that constitutes our model for the S-shaped component. The symmetry axis of the 
precessing jet is indicated by the thick dashed line, which forms an angle i to the plane of the 
sky and lies at a position angle ip in the plane of the sky. The amplitude of the precessing jet, 
al, increases with the distance / away from the driving source. The jet precesses on a cone 
with a half-opening angle of 6, which is represented by two dash-dot lines. The red and blue 
cones represent the outflow cavity traced by the X-shaped component (see Fig. 8), which 
has the opposite velocity pattern to the S-shaped component. The blue/red arrows in the 
upper-right insert shows the directions in which a bent jet emerges from its driving source 
in the model discussed in §4.3.2. The black dashed line is orthogonal to the circumstellar 
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Fig. 11. — Predicted morphology (upper panels) and line-of-siglit velocity along the path 
(lower panel) of a precessing outflow, driven by a straight jet, with model parameters listed 
in the first row of Table 4. The black curves in the upper panels show the path of the jet 
superposed on the integrated CO (2-1) intensity map of Figure 4 with the central compact 
component included (left panel) and with it subtracted out (right panel). The black curves 
in the lower panels show the derived line-of-sight velocity of the outflow along the path of the 
jet superposed on the central C0(2-l) velocities derived from gaussian fits to the spectra at 
regularly spaced intervals along this path (see text). The smaller solid error bars correspond 

tn mfiasiirpmpnt imr.firtaintifis in thfi rp.ntral vp.lnnitifis and thp, larTOr larTO da.c;hfid prrnr bars 
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Fig. 12. — Same as for Figure 11, but in this case for a precessing outflow driven by a bent jet 
with the model parameters listed in the second row of Table 4. This model better reproduces 
the observed line-of-sight velocities as shown in the lower panel. 
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Fig. 13. — Same as for Figure 11, but in this case for an outflow driven by a mirror-symmetric 
jet with model parameters as described in §4.3.3. 



